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In-vitro  retinal  model 
reveals  a  sharp  transition 
between  laser  damage 
mechanisms 

Michael  L.  Denton,3  C.  D.  Clark  III,3  Michael  S.  Foltz,3 
Kurt  J.  Schuster,3  Gary  D.  Noojin,3  Larry  E.  Estlack,b 
and  Robert  J.  Thomasc 
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bConceptual  MindWorks,  Inc.,  San  Antonio,  Texas  78228 
cAir  Force  Research  Laboratory,  711  HPW/RHDO, 

Brooks  City-Base,  Texas  78235 

Abstract.  We  use  laser  damage  thresholds  in  an  in-vitro 
retinal  model,  and  computational  simulations  to  examine 
the  laser  exposure  durations  at  which  damage  transitions 
from  photothermal  to  photochemical  at  413  nm.  Our  re¬ 
sults  indicate  a  dramatic  shift  in  1 -h  damage  thresholds 
between  exposure  durations  of  60  and  1 00  s.  The  trend  in 
our  in-vitro  results  is  similar  to  a  trend  found  in  a  recent 
study  where  retinal  lesions  were  assessed  1 -h  post  laser 
exposure  in  the  rhesus  eye  Our  data  suggest  that  nonther- 
mal  mechanisms  did  not  significantly  contribute  to  cell 
death,  even  for  exposures  of  60  s.  Knowledge  of  the  tran¬ 
sition  point,  and  lack  of  concurrent  thermal  and  nonther- 
mal  damage  processes,  are  significant  for  those  wishing  to 
devise  a  comprehensive  computational  damage  model. 
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Cells  exposed  to  laser  pulses  longer  than  the  thermal  dif¬ 
fusion  rate  can  experience  damage  by  either  photothermal  or 
photochemical  mechanisms,  depending  on  wavelength,  expo¬ 
sure  duration,  and  beam  diameter.  In  addition,  laser  irradiance 
(E)  dictates  the  degree  of  chromophore-dependent  heating  in 
tissues.  In  general,  the  irradiance  requirement  for  damage  is 
reduced  as  exposure  duration  is  extended.  This  trend  explains 
why  wavelength-dependent  nonthermal  damage  only  pre¬ 
dominates  for  long  laser  exposures,  and  how  thermal  damage 
can  be  forced  to  occur  from  relatively  short  exposures.  Radi¬ 
ant  exposure  (H)  is  used  to  correlate  laser  dose  to  damage 
outcomes,  and  is  calculated  by  multiplying  irradiance 
(W  cm-2)  by  exposure  duration. 

In  animals,  where  measurements  of  temperature  rise  and 
oxidative  products  are  difficult,  researchers  have  relied  on 
evaluations  of  trends  in  damage  thresholds  relative  to  wave¬ 
length  (action  spectrum)  and  exposure  duration  [temporal  ac¬ 
tion  profile  (TAP)]1  to  distinguish  a  damage  mechanism.  The 
principal  of  irradiance  reciprocity,  where  threshold  irradiance 


Address  all  correspondence  to:  Robert  J.  Thomas,  E-mail: 

robert.thomas@brooks.ai.mil 

Journal  of  Biomedical  Optics 


is  inversely  proportional  to  exposure  duration,  is  a  prominent 
feature  of  photochemical  (nonthermal)  damage,  and  is  easily 
identified  in  TAP  analyses  using  threshold  irradiance  (ETAP) 
or  threshold  radiant  exposure  (HTAP).  Ocular  lesions  as¬ 
sessed  at  24-h  post  exposure  (or  longer)  for  exposures  of 
100  s  (or  longer)  at  514  nm  (or  shorter)  were  considered 
photochemical.1'2  When  lesions  from  514-nm  exposures  were 
assessed  at  1  h,  damage  from  a  1000-s  exposure  was  consis¬ 
tent  (HTAP)  with  photothermal  processes.2  Therefore,  this 
“latency”  in  the  appearance  of  damage  is  also  useful  in  char¬ 
acterizing  laser  damage  mechanisms. 

While  characterizing  a  novel  in-vitro  retinal  model,3  we 
made  note  of  a  peculiar  trend  for  damage  thresholds  at 
413  nm  that  compared  with  in-vivo  thresholds  at  457.9  and 
441.6  nm  reported  by  Lund,  Stuck,  and  Edsall.4  An  HTAP 
analysis  of  the  data  revealed  thresholds  that  continued  along 
the  thermal  trend  line  to  threshold  radiant  exposure  values 
greater  than  that  for  100  s,  which  was  expected  to  be  at  irra¬ 
diance  reciprocity.  Figure  1(a)  is  an  HTAP  for  the  441.6-nm 
(327-/xm  beam  diameter)  data  taken  from  Ref.  4.  Without  the 
usual  log-log  axis,  the  HTAP  clearly  shows  that  the  1  -h  post 
exposure  threshold  values  for  the  5-,  8-,  and  16-s  exposures 
are  greater  than  that  of  the  photochemical  threshold 
(26  J  cm-2).  As  the  authors  concluded,4  a  lack  of  latent  le¬ 
sions  following  laser  exposures  of  5  s  and  shorter  indicated 
damage  by  thermal  mechanisms.  However,  the  fact  that 
threshold  values  for  the  8-  to  16-s  exposures  did  not  decrease 
to  the  level  of  the  100-s  threshold  implied  the  possibility  of 
mixed  damage  mechanisms  in  the  progression  to  purely  pho¬ 
tochemical  processes.  The  authors  discussed  how  the  16-s 
threshold  data  suggested  a  gradual  transition  between  mecha¬ 
nisms  due  to  a  delayed  appearance  of  threshold  lesions  (indi¬ 
cating  nonthermal  damage)  and  a  dependence  on  the  diameter 
of  the  irradiated  retinal  area  (indicating  thermal  damage). 

In  this  study,  we  use  a  combination  of  new  in-vitro  damage 
data  and  computational  simulations  of  both  temperature  rise 
and  damage  thresholds  at  413  nm  to  address  the  shift  to  pho¬ 
tochemical  damage  processes  in  more  detail.  Except  for  the 
following  changes,  laser  exposures  were  as  described 
previously.'1  Retinal  pigment  epithelial  (RPE)  cells  (about  160 
melanosomes/cell)  were  exposed  to  413-nm  laser  irradiation 
in  48-well  plates  containing  0.1-mL  complete  Hank’s  bal¬ 
anced  salt  solution.  Cells  were  exposed  to  a  0.3-mm-diam 
flat-top  beam  (via  an  88-mm  focal  length  lens)  in  an  environ¬ 
mentally  controlled  enclosure  (Fig.  2)  that  ensured  consistent 
temperature  (35  to  37  °C)  and  relative  humidity  (60  to  70%). 
Estimated  dose  for  50%  lethality  (ED50)  values  were  calcu¬ 
lated  using  the  Probit  method,5'6  where  Probit  slopes  represent 
the  first  derivative  with  respect  to  dose  at  a  probability  of  0.5. 
Systematic  uncertainty  (15%)  in  our  irradiance  values  was 
calculated  as  previously  described.7 

To  simulate  the  experimental  exposures,  we  used  a  laser- 
tissue  damage  program  developed  by  our  group.8  The  model 
uses  a  laser  propagation  model  (geometric  optics  assumed)  to 
compute  a  source  term  (implementing  Beer’s  law  of  linear 
absorption)  for  a  thermal  heat  solver.  The  heat  solver  accounts 
for  the  multiple  layers  comprising  the  system  (buffer,  cells, 
and  plastic  well  plate)  and  air/surface  boundary  conditions  to 
predict  temperature  rise  at  the  7-/rm-thick  cell  layer.  These 
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Fig.  1  Comparison  of  in-vivo  and  in-vitro  damage  trends,  (a)  HTAP  for  1  -  and  48-h  damage  thresholds  (441 .6  nm)  in  the  in-vivo  MVL  study  of  Ref. 
4.  The  100-s  48-h  data  point  is  expected  to  represent  irradiance  reciprocity,  (b)  HTAP  for  1-h  damage  thresholds  (413  nm)  in  the  in-vitro  retinal 
model.  Estimated  thresholds  using  the  BTEC  thermal  model  are  represented  by  the  dashed  line.  Error  bars  depict  upper  and  lower  fiducial  limits  (95 
%  confidence  intervals). 


temperatures  were  input  into  an  Arrhenius  rate  equation  that 
was  numerically  integrated  to  determine  damage  outcomes. 
The  two  Arrhenius  rate  parameters  (A  =  3.1  X  10"  s_1;  Ea 
=6.28  X  105  J/mol_1)  used  in  the  damage  integral  (11)  calcu¬ 
lation  were  those  reported  by  Welch  and  Polhamus.9  A  nu¬ 
merical  search  algorithm  was  used  to  determine  the  threshold 
irradiance  that  solved  for  an  fl  value  of  1  at  the  center  of  the 
beam. 

The  1-h  ED50  irradiances  for  the  current  in-vitro  exposures 
are  given  in  Table  1.  Consistent  environmental  conditions  dur¬ 
ing  exposures  led  to  low  variance  (fiducial  limits)  about  the 
threshold  values.  Notice  that,  as  expected,  the  irradiance  re¬ 
quirement  for  damage  was  reduced  by  extending  the  duration 
of  the  laser  exposure,  although  there  was  no  significant 
change  between  the  40-  and  60-s  thresholds.  Table  1  also 
shows  that  the  threshold  value  for  the  200-s  exposure  was 
exactly  half  the  100-s  threshold,  defining  irradiance  reciproc¬ 
ity  (nonthermal  damage)  for  the  0.3-mm  beam  at  36  °C  in 
the  in-vitro  retinal  system. 

An  HTAP  analysis  of  the  in-vitro  results  [Fig.  1(b)]  shows 
a  trend  similar  to  the  1-h  assessments  of  Ref.  4  at  441.6  nm 
[Fig.  1(a)].  Power  functions  describing  the  in-vitro  and  in-vivo 
data  had  nearly  identical  exponents  (0.76  and  0.78,  respec¬ 


Fig.  2  Schematic  representation  of  41  3-nm  beam  delivery  for  in-vitro 
laser  exposures.  M,  mirror;  X/2,  half-wave  plate;  Pol,  polarizing 
beamsplitter;  T,  optical  telescope;  BSh,  beam  shaper;  BS,  beamsplitter; 
S,  shutter;  L,  lens;  ND,  neutral  density  filter;  CCD,  charge-coupled 
device  camera;  and  Obj,  microscope  objective. 


tively),  but  the  in-vivo  curve  was  shifted  to  lower  threshold 
values  by  a  factor  of  about  6.5,  presumably  due  to  differences 
in  melanosome  density  as  discussed  previously.3  Although 
difficult  to  see  in  Fig.  1(b),  the  threshold  data  points  for  the 
0.1-  and  1.0-s  exposures  are  near  the  origin.  Notice  the  in¬ 
verse  relationship  between  irradiance  and  radiant  exposure, 
such  that  the  0.1  -s  exposure  required  the  greatest  irradiance 
(Table  1)  and  the  least  radiant  exposure  [Fig.  1(b)]  for  gener¬ 
ating  threshold  damage.  However,  as  seen  in  Fig.  1,  this  trend 
is  dramatically  broken  when  the  damage  mechanism  shifted  to 
nonthermal,  as  indicated  by  irradiance  reciprocity.  This  nec¬ 
essarily  means  that  there  was  a  significant  thermal  component 
in  the  death  process  for  exposures  of  1  min.  The  correlation 
coefficient  for  the  power  function  describing  the  0.1-  to  60-s 
data  [see  Fig.  1(b)]  would  suggest  a  similar  degree  of  thermal 
component  over  this  entire  range  of  exposure  durations. 

As  a  means  of  understanding  the  thermal  component  of  the 
in-vitro  damage  thresholds,  we  simulated  irradiance  threshold 
values  (f!=l)  for  exposure  durations  of  0.1  to  200  s  and 

Table  1  1-h  damage  threshold  irradiance  values  for  41 3-nm  laser 
exposures  of  in  in-vitro  retinal  model.  FL  are  Fiducial  limits  (95% 
confidence  interval).  Systematic  uncertainty  was  determined  to  be 
1  5%. 


Exposure 
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0.1 

94 
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92 
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18 
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41.3 
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18 

40 
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38.0 

15 
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9 
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Fig.  3  S  imulated  (BTEC  model)  peak  temperature  rise  values  for  in- 
vitro  ED50  irradiance  values  (circles)  as  a  function  of  exposure  dura¬ 
tion.  Peak  temperatures  of  irradiances  extrapolated  in  accordance 
with  reciprocity  (t-1)  were  also  simulated  and  shown. 


plotted  the  radiant  exposure  results  alongside  the  in-vitro  data 
[Fig.  1(b)].  The  simulated  curve  fits  fairly  well  with  the  in- 
vitro  data,  falling  within  the  95%  confidence  intervals  for  all 
but  the  40-s  data  point.  The  deviation  in  simulated  and  in- 
vitro  HTAP  curves  suggests  that  mixed  damage  mechanisms 
may  exist  for  exposures  of  40  to  60  s.  However,  we  under¬ 
stand  that  our  choice  in  values  for  the  Arrhenius  rate  param¬ 
eters,  the  fact  that  one  of  these  parameters  (frequency  factor 
A)  has  a  slight  dependence  on  temperature,  and  our  15%  sys¬ 
tematic  uncertainty  (Table  1),  are  all  factors  that  could  lead  to 
the  observed  disparity. 

Looking  at  the  irradiance  thresholds  in  Table  1,  we  see  no 
difference  between  the  40-  and  60-s  exposures,  which  means 
the  continuation  of  radiant  exposure  thresholds  along  the 
power  curve  in  the  HTAP  is  directly  proportional  to  the  in¬ 
crease  in  exposure  duration.  However,  this  is  not  necessarily 
unexpected,  because  threshold  irradiance  (and  temperature)  is 
proportional  to  the  inverse  of  log  time  under  the  Arrhenius 
damage  model,  which  predicts  a  convergence  of  threshold  ir¬ 
radiances  (and  temperatures)  at  longer  exposure  durations.  A 
similar  scenario  appears  to  play  out  in  the  rhesus  eye.  The 
tabulated  retinal  1-h  threshold  irradiance  values  for  the 
5-  to  16-s  exposures  reported  in  Ref.  4  were  all  three  statis¬ 
tically  identical  as  well. 

Previous  authors  have  used  a  peak  temperature  rise  of 
9  to  10  °C  to  signify  a  minimum  requirement  for  damage  by 
thermal  means.10  As  expected,  simulated  peak  temperature 
rises  (Fig.  3)  for  the  100-  and  200-s  in-vitro  threshold  values 
were  low  (3.8  and  1.9  °C,  respectively).  Figure  3  also  shows 
that  simulated  peak  temperature  rises  of  greater  than  10  °C 
were  calculated  for  the  in-vitro  ED50  irradiance  values  for 
exposures  of  60  s  and  shorter,  which  implies  that  sufficient 
heat  was  generated  to  produce  death  by  thermal  mechanisms. 
When  we  looked  at  the  individual  damage  outcome  data 
(damage  versus  no  damage)  for  the  60-s  exposures,  we  found 
that  the  lowest  irradiance  (23.3  W  cm-2)  that  caused  damage 
corresponded  to  a  (simulated)  temperature  rise  of  9.3  °C. 
Again,  this  indicated  that  all  the  60-S  exposure  damage  out¬ 
comes  had  sufficient  temperature  rises  to  cause  damage  by 
thermal  means.  However,  this  does  not  exclude  the  possibility 
of  nonthermal  (photo-oxidation)  events  occurring  concur¬ 
rently  with  this  elevation  of  temperature.  The  lack  of  an  in¬ 
termediate  temperature  rise  for  the  40-  and  60-s  ED50  irradi¬ 


ance  values,  such  as  5  to  9  °C,  suggests  that  if  concurrent 
mechanisms  did  exist,  they  were  neither  additive  nor  syner¬ 
gistic  to  the  overall  damage  rate  process. 

Finally,  we  were  interested  in  the  expected  temperature 
rise  of  cells  receiving  the  40-  and  60-s  reciprocity  irradiance 
doses.  On  extrapolating  it~v)  from  the  100-  and  200-s  irradi¬ 
ances,  we  simulated  peak  temperatures  of  9.4  and  6.3  °C  for 
the  extrapolated  irradiances  corresponding  to  40-  and  60-s 
exposures,  respectively.  This  suggests  that  if  there  did  exist 
additive  or  synergistic  effects  from  thermal  and  nonthermal 
damage  mechanisms,  it  would  be  manifested  in  the  data  for 
60-s  exposures. 

In  conclusion,  the  in-vitro  retinal  model  showed  transitions 
in  the  damage  mechanism  for  1-h  thresholds  similar  to  those 
found  previously  in  an  in-vivo  model.  In  the  in-vitro  model, 
the  transition  from  photothermal  to  nonthermal  damage  was 
sudden,  occurring  somewhere  between  exposure  durations  of 
60  and  100  s.  Additional  data  in  this  exposure  range  are  being 
collected  to  characterize  this  transition.  The  current  data  can¬ 
not  rule  out  the  possibility  of  photochemical  oxidation  occur¬ 
ring  during  damaging  60-s  exposures  at  413  nm.  We  are  cur¬ 
rently  conducting  experiments  to  address  this  issue.  However, 
we  believe  that  nonthermal  processes  do  not  contribute  to  cell 
death  at  1-h  postexposure. 
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